Spin Waves and Electronic Interactions in La2Cu04 
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The magnetic excitations of the square-lattice spin-1/2 antiferromagnet and high-Tc parent 
La2CuC>4 are determined using high-resolution inelastic neutron scattering. Sharp spin waves with 
absolute intensities in agreement with theory including quantum corrections are found throughout 
the Brillouin zone. The observed dispersion relation shows evidence for substantial interactions 
beyond the nearest-neighbor Heisenberg term, which can be understood in terms of a cyclic or ring 
exchange due to the strong hybridization path around the CU4O4 square plaquettes. 
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While there is consensus about the basic phenomenol- 
ogy - electron pairs with non-zero angular momentum, 
unconventional metallic behavior in the normal state, 
tendencies towards inhomogeneous charge and spin den- 
sity order - of the high temperature copper oxide super- 
conductors, there is no agreement about the microscopic 
mechanism. After over a decade of intense activity, there 
is not even consensus as to the simplest "effective Hamil- 
tonian" , which is a short-hand description of the mo- 
tions and interactions of the valence electrons, needed 
to account for cuprate superconductivity. Because much 
speculation is centered on magnetic mechanisms for the 
superconductivity, it is important to identify the inter- 
actions among the spins derived from the unfilled Cu 2+ 
d-shells. The present experiments show that there are 
significant (on the scale of the pairing energies for high- 
Tc superconductivity) interactions coupling spins at dis- 
tances beyond the 3.8 A separation of nearest-neighbor 
Cu 2+ ions. Cyclic or ring exchange due to a strong hy- 
bridization path around the CU4O4 squares (see Fig. |]A), 
from which the cuprates are built, provides a natural ex- 
planation for the measured dispersion relation. Cu02 
planes are thus the second example of an important Fermi 
system ( 3 He is the other Q) where significant cyclic ex- 
change terms have been deduced. 

Magnetic interactions are revealed through the 
wavevector dependence or dispersion of the magnetic ex- 
citations. In magnetically ordered materials, the domi- 
nant excitations are spin waves which are coherent (from 
site to site as well as in time) precessions of the spins 
about their mean values. The lower frame of Fig. |l|B 
shows the dispersion relation calculated using conven- 
tional linear spin-wave theory in the classical large-.!? 
limit, where the only magnetic interaction is a strong 
nearest- neighbor superexchange coupling J Q . We iden- 
tify wavevectors by their coordinates (h, k) in the two- 
dimensional (2D) reciprocal space of the square lattice. 



Spin waves emerge from the wavevector (1/2,1/2) charac- 
terizing the simple antiferromagnetic (AF) unit cell dou- 
bling in La2Cu04 ||, and disperse to reach a maximum 
energy 2 J that is a constant along the AF zone bound- 
ary marked by dashed squares in Fig. [l]B. Longer-range 
interactions manifest themselves most simply at the zone 
boundary. The upper frame of Fig. |l|B shows the disper- 
sion calculated with modest interactions between next 
nearest- neighbors. Virtually the only visible effect of the 
additional interactions is the dispersion of the spin waves 
along the zone edge. Thus, experiments to test for such 
interactions must measure the spin waves along the zone 
boundary. Only inelastic neutron scattering with high 
energy and wavevector resolution can accomplish this, 
although photon spectroscopy j|-f7| has led to suspicions 
of such interactions. 

For La2CuC>4, a requirement that complicates meeting 
the resolution goals is the need to use neutrons with en- 
ergies in the epithermal, 0.1-1.0 eV, range rather than in 
the more conventional cold and thermal 2-50 meV, 
regimes. An early high energy neutron scattering ex- 
periment H revealed well-defined spin-wave excitations 
throughout the Brillouin zone which could be modeled 
using a nearest-neighbor Heisenberg exchange J=136 
meV. The directions of the scattered neutrons were spec- 
ified only to within the solid angle determined by the 
large detector dimensions. Thus, the measured spectra 
represented averages over large portions of the reciprocal 
space, so that dispersion along the zone boundary was 
unresolvable and only an upper bound could be placed 
on further neighbor couplings. The advance enabling the 
present investigation is the use of position-sensitive de- 
tectors for the scattered neutrons, which increases the 
wavevector resolution by an order of magnitude. The new 
detector bank is installed in the direct-geometry High- 
Energy Transfer (HET) time-of-flight spectrometer at the 
ISIS proton-driven pulsed neutron spallation source. 
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FIG. 1. (color) A The CuCb plane showing the atomic 
orbitals (Cu 3d x 2_ y 2 and O 2p x ,y) involved in the mag- 
netic interactions. J, J' and J" are the first-, second- and 
third-nearest-neighbor exchanges and J c is the cyclic interac- 
tion which couples spins at the corners of a square plaquette. 
Arrows indicate the spins of the valence electrons involved in 
the exchange. B Lower surface is the dispersion relation for 
J— 136 meV and no higher-order magnetic couplings or quan- 
tum corrections. The upper surface shows the effect of the 
higher-order magnetic interactions determined by the present 
experiment. Color is spin-wave intensity. 

Fig. |2|A shows data in the form of constant energy 
scans for wavevectors around the antiferromagnetic zone 
center. As E increases, counter-propagating modes be- 
come apparent. As the zone boundary is approached and 
there is less dispersion, inspection of Fig. ^]B reveals that 
it should be easier to locate the spin waves via energy 
scans performed at fixed wavevector. Fig. shows a 
series of such scans collected at various points along the 
zone boundary. The spin waves have a clearly notice- 
able dispersion, from a minimum of 292±7 meV near 
Q=(3/4,l/4) to a maximum of 314±7 meV near (1/2,0). 
This is in obvious contrast to the dispersion-less behav- 
ior of linear spin-wave theory for the nearest-neighbor 
Heisenberg model. We have collected data throughout 
the Brillouin zone and Fig. ||A shows the resulting dis- 
persion along major symmetry directions obtained from 
cuts of the type shown in Fig. ||. Fig. |^B displays the cor- 
responding spin- wave intensities, in absolute units cali- 
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FIG. 2. Scattering from the spin waves in La2Cu04 
(T=295 K). Data result from 68 hours (A) and 98 hours (B) 
counting at a proton current of 170 /j,A. The sample is de- 
scribed in Jl^] . Solid lines are fits to a spin- wave cross-section 
convolved with the instrumental resolution. (A) Const-iJ cuts 
near the AF zone center for an incident energy Si=250 meV. 
Q z wavevector components at scan centers are Z=2.8 (bottom 
panel), 6.2 and 9.6 r.l.u. Open circles are a background mea- 
sured near the (0,0) position. Dashed curve is the instrumen- 
tal response to spin waves of infinite velocity. (B) Const-Q 
cuts, with Ei=750 meV, yield the dispersion along the AF 
zone boundary. Vertical dotted line at £=300 meV is a guide 
to the eye. I values at peak position vary from 8.8 (bottom 
panel) to 9.5 (top panel). A background measured near the 
nuclear zone center (1,0) has been subtracted. Dashed curve 
is the instrumental response to a dispersionless mode. 



brated using acoustic phonon scattering from the sample. 

To understand our results, we consider a Heisenberg 
Hamiltonian including higher order couplings [Rl] O] 



+j c ]T {(s i -s j )(s k -Si) + (s i -s l )(s k -s j ) 

(i,j,k,l) 

-(Si-SkXSj-Si)}, (i) 

where J, J' and J" are the first-, second- and third- 
nearest-neighbor magnetic exchanges where the paths 
are illustrated in Fig. ^A. J c is the ring exchange in- 
teraction coupling four spins (labelled clockwise) at the 
corners of a square plaquette. Each spin coupling is 
counted once in Eq. |l]). Using classical (large-5) lin- 
ear spin- wave theory the dispersion relation is [jl3jl5 

uj q =2Z c {Q)^A 2 q -B 2 Q , A Q =J - J c /2 - (J' - J c /4)(1 - 

u h v h )-J" [1 - {v 2h + v 2k )/2], B Q ={J-J c /2){v h +v k )/2, 
v x = cos(27rx) and Z C (Q) is a renormalization factor |l6| ] 
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that includes the effect of quantum fluctuations. Within 
linear spin-wave theory all three higher-order spin cou- 
plings (J 1 , J" and J c ) have similar effects on the dis- 
persion relation and intensity dependence, therefore they 
cannot be determined independently from the data with- 
out additional constraints. We first assume that only 
J and J' are significant as in [Jf7|], i.e. J"=J c =0. The 
solid lines in Fig. ^ are fits to a one-magnon cross-section 
and Fig. || shows fits to the extracted dispersion rela- 
tion and spin-wave intensity. As can be seen in the fig- 
ures, the model provides an excellent description of both 
the spin-wave energies and intensities. The extracted 
nearest- neighbor exchange J=111.8±4 meV is antiferro- 
magnetic, while the next-nearest-neighbor exchange J'=- 
lf.4±3 meV across the diagonal is ferromagnetic. A 
wavevector-independent quantum renormalization factor 
|l6f Z c =1.18 was used in converting spin- wave energies 
into exchange couplings. The zone-boundary dispersion 
becomes more pronounced upon cooling as shown in Fig. 
HA and the dispersion at T—10 K can be described by 
the couplings J=104.1±4 meV and J'=-18±3 meV. 

A ferromagnetic J' contradicts theoretical predictions 
p9| , which give an antiferromagnetic superexchange J'. 
Wavevector-dependent quantum corrections |20| to the 
spin-wave energies can also lead to a dispersion along 
the zone boundary even if J' = 0, but with sign oppo- 
site to our result. Another problem with a ferromagnetic 
J' comes from measurements on Sr2Cu304Cl2 This 
material contains a similar exchange path between Cu 2+ 
ions to that corresponding to J' in La2CuC>4 and analysis 
of the measured spin- wave dispersion leads to an antifer- 
romagnetic exchange coupling for this path [ pl| . 

While we can cannot definitively rule out a ferromag- 
netic J' we can obtain a natural descri pti on of the data 
in terms of a one-band Hubbard model P2| , an expansion 
of which yields the spin Hamiltonian in Eq. ([!]) where the 
higher-order exchange terms arise from the coherent mo- 
tion of electrons beyond nearest-neighbor sites jn] Q| . 
The Hubbard Hamiltonian has been widely used as a 
starting point for theories of the cuprates and is given by 
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H.c. 



(2) 



where stands for pairs of nearest-neighbors counted 
once. Eq. (||) has two contributions: the first is the 
kinetic term characterized by a hopping energy t be- 
tween nearest-neighbor Cu sites and the second the po- 
tential energy term with U being the penalty for dou- 
ble occupancy on a given site. At half filling, the case 
for La2CuC>4, there is one electron per site and for 
t/U — > 0, charge fluctuations are entirely suppressed 
in the ground state. The remaining degrees of freedom 
are the spins of the electrons localized at each site. For 
small but non-zero t/U, the spins interact via a series of 
exchange terms, as in Eq. (Jl|) , due to coherent electron 
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FIG. 3. (A) Dispersion relation along high symmetry di- 
rections in the 2D Brillouin zone, see inset (C), at T=10 K 
(open symbols) and 295 K (solid symbols). Squares were ob- 
tained for Ei=250 meV, circles for i?i=600 meV and triangles 
for Ei=750 meV. Points extracted from constant-iJ(-Q) cuts 
have a vertical (horizontal) bar to indicate the E(Q) inte- 
gration band. Solid (dashed) line is a fit to the spin-wave 
dispersion relation at T—10 K (295 K) as discussed in the 
text. (B) Wavevector-dependence of the spin-wave intensity 
at T=295 K compared with predictions of linear spin-wave 
theory shown by the solid line. The absolute intensities 
yield a wavevector-independent intensity-lowering renormal- 
ization factor of 0.51±0.13 in agreement with the theoretical 
prediction of 0.61 jh| that includes the effects of quantum 
fluctuations. 



motion touching progressively larger numbers of sites. 
If the perturbation series is expanded to order t A (i.e. 
4 hops) one regains the Hamiltonian ([!]) with the ex- 
change constants J = At 2 /U - 24t i /U 3 , J c = 80i 4 /C/ 3 
and J' = J" = At 4 /U 3 jn]-||. We again fitted the 
dispersion and intensities of the spin-wave excitations 
using these expressions for the exchange constants and 
linear spin wave theory. The fits are indistinguishable 
from those for variable J and J'. Again assuming |2^] 
Z c =1.18, we obtained t=0.33±0.02 eV and C/=2.9±0.4 
eV (T—295 K), in agreement with t and U determined 
from photoemission |24|] and optical spectroscopy p5[ . 
The corresponding exchange values are J = 138.3±4 meV, 
J c =38±8 meV and J'=J"= J c /20=2±0.5 meV (the pa- 
rameters at T=10 K are <=0.30±0.02 eV, C/=2.2±0.4 eV, 
J=146.3±4 meV and J c =61±8 meV). Using these values, 
the higher-order interactions amount to ~11% (T=295 
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K) of the total magnetic energy 2(J-J c /4-J'-J") required 
to reverse one spin on a fully-aligned Neel phase. 

Many results on oxides of copper fall into place when 
cyclic exchange of the size extracted from our experi- 
ments is taken into account. First, the relative magni- 
tude of the cyclic exchange J c / J=0.27±0.06 at T=295 
K (0.41±0.07 at T=10 K) is similar to the ratio of 



0.30 estimated from numerical simulations |26| on finite 
clusters taken from the Cu-0 square lattice. Second, 
magnetic Raman scattering [^|j5j and infrared absorption 
experiments ||[7]] show an unusual broadening towards 
higher energies that cannot be accounted for by a sim- 
ple (quadratic) Heisenberg Hamiltonian, but can be at- 
tributed d|| to a cyclic term. Finally, in the related com- 
pound Sri4Cu2404i, which has square plaquettes stacked 
to form a ladder, the exchange constants corresponding 
to the nearly-equal-length rungs and legs of the ladder 
are 130 meV and 72 meV respectively p]J when no cyclic 
exchange is included. The inclusion of a ring exchange 
term [H J c =34 meV allows the rungs and legs of the 
ladder to have similar exchange constants of 121 meV. 

We have used a new high-wavevector-resolution 
epithermal-neutron scattering technique to discover that 
interactions beyond those coupling nearest-neighbor 
Cu 2+ ions are needed to account for the magnetism 
of La2CuC>4. The observed further neighbor couplings 
may be explained by a four-spin cyclic interaction, which 
arises because the large orbital hybridization in the Cu02 
planes provides an exchange path to include all four spins 
at the corners of elementary CU4O4 plaquettes. Thus, 
La2CuC>4 joins the nuclear magnet 3 He pj as a sys- 
tem where there is good evidence for substantial ring 
exchange. Ring exchange occurs even for very simple 
models, such as the single band Hubbard model, which 
contains only hopping (t) and on-site Coulomb (U) terms. 
We have determined t and U using only knowledge - in 
the form of our spin wave dispersion relation - of charge- 
neutral excitations and find values in excellent agreement 
with those obtained via charge-sensitive spectroscopies as 
well as numerical work on finite clusters. Thus, our re- 
sults demonstrate that a one-band Hubbard model is an 
excellent starting point |2j| for describing the magnetic 
interactions in the cuprates, and that even when consid- 
ering relatively low energy spin excitations, charge fluc- 
tuations involving double occupancy must be taken into 
account. In addition, the scale of the cyclic exchange in- 
teractions, which are comparable to pairing energies in 
the high-Tc materials, implies that they themselves or 
related electronic currents po|| might be important for 
superconductivity in the doped cuprates. 
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